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The contribution of the proton polarizability to the ground state hyperfine splitting in muonic 
hydrogen is evaluated on the basis of modern experimental and theoretical results on the proton 
polarized structure functions. The value of this correction is equal to 4.6(8) • 10 -4 times the Fermi 
splitting. 

36.10 Dr, 12.20 Ds, 31.30 Jv 

Keywords: muonic hydrogen, hyperfine splitting, proton polarizability, nucleon polarized structure functions 



I. INTRODUCTION 



The investigation of the energy spectra of hydrogenic atoms (positronium, muonium, hydrogen atom, muonic 
hydrogen, et. al) as well as lepton anomalous magnetic moments provides a test of quantum electrodynamics and the 
theory of electromagnetic bound states with very high accuracy. Moreover, the values of the fundamental physical 
constants (the particle masses, fine structure constant, proton charge radius, etc.) can be determined more precisely. 
Inclusion of new simple atomic systems in the range of the experimental investigations can lead to the significant 
progress in solving of these problems. The measurement of the muonic hydrogen Lamb shift at PSI with a precision 
of 30 ppm will allow to improve our knowledge of the proton charge radius by an order of magnitude Q . Another 
important problem is connected with the measurement of the ground state hyperfine splitting (HFS) in muonic 
hydrogen || . In the case of electronic hydrogen HFS was measured with extremely high accuracy many years ago [g| : 

A ^HFs( e P) = 1420405.7517667(9) kHz. (1) 

The corresponding theoretical expression of the hydrogen hyperfine splitting can be written in the form 
(AE^s = 27rhA^ h FS ) §: 

AE&s = E F (1 + 5^ + S S + O, E F = | a 4 Mpmjml 

3 (mi + m 2 )' i 

where /ip is the proton magnetic moment, mi, 1112 are the masses of the electron (muon) and proton. The calculation of 
different corrections to Ep has a long history. Modern status in the theory of hydrogenic atoms was presented in details 
in (5Q ED denotes the contribution of higher-order quantumelectrodynamical effects. This correction is known with 
an accuracy 10~ 7 Q. Corrections 8 s and S p take into account the influence of strong interaction. 8 s describes the 
effects of the proton finite-size and recoil contribution. S p is the correction due to the proton polarizability. The 
calculation of the nucleon electromagnetic form factors and the polarized structure functions on the basis of QCD 
inspired effective field theories || with high degree of accuracy can lead to the determination of the contributions 
8 and 8 F in the hydrogen HFS. The progress in solving of these problems strongly depends also on the successful 
realization of the program for the study of the nucleon excitations at CEB AF || . 

The comparison of the theoretical result (2) for electronic hydrogen without the unknown proton polarizability 
correction with the experimental value (1) yields: 
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where the main sources of uncertainty in (3) are the inaccuracy of the proton form factor parameterization (dipole fit 
etc.) and the contradictory experimental data on the proton radius. The general structure of the expression for the 
muonic hydrogen hyperfine splitting coincides with (2). The QED contribution emerges after replacement m — > m^. 
The corrections 8 (/^p) and 5 p (/xp) also depend on the lepton mass. So, there are no simple relations between 8 and 
<5 P for muonic and electronic hydrogen. But the inclusion of the muonic hydrogen HFS in the sphere of the theoretical 
and experimental investigations can shed additional light on the proton polarizability and proton structure corrections 
in both systems. 

II. GENERAL FORMALISM 

The main part of the one-loop proton structure correction is determined by the following expression (Zemach 
correction) @,|]: 



ae = e f 2 ^ r dl5 



7T 2 J ( P 2 +b 2 ) 2 



Ge(-p 2 )G m (-P 2 ) 1 

1 + K 



= E F (-2^a)Rp, b = a^, (4) 



where R p is the Zemach radius, /i is reduced mass. In the coordinate representation the Zemach correction is 
determined by the magnetic moment density Pm( t ) and by the density of electric charge Pe( t )- The value of R p can 
be obtained in the analytical form by using the dipole parameterization for the proton electromagnetic form factors 
Ge and Gm: 



35 _ 126 1056 2 fb 3 \ 
8A ~ A 2 " + ~4A r + V A 1 J ' 



R P = — - — + — + o{-), (5) 



where A = 0.898m p . Correction to the mass independent part 35/8A in eq. (5) of order 0(b/A) is small due to the 
first power of the fine structure constant. Relative order contribution of the Zemach correction to the hydrogen HFS 
is the following: 

electronic hydrogen : R p = 1.025 fm, <5(Zemach) = —38.72 ppm, (6) 

muonic hydrogen : R p = 1.022 fm, <5(Zemach) = -71.80 • 10~ 4 , (7) 

Other possible parameterizations for the nucleon electromagnetic form factors Q lead to changing the Zemach correc- 
tion by 3-4 %. The second recoil part of the one-loop contribution to the HFS in hydrogen is equal to 5.22 • 10~ 6 E F |§. 
Other uncertainty of the expression (2) is connected with the proton polarizability correction 8 P . This contribution 
can be obtained from the two-photon exchange diagrams as illustrated in Fig. 1. The corresponding amplitudes of the 
virtual Compton scattering on the proton can be expressed through nucleon polarized structure functions G\(v, Q~ ) 
and G 2 (z^, Q 2 ). The inelastic contribution of the diagrams (a), (b) Fig. 1 to the hydrogen HFS was studied in [[To] pT[ , 
where the electron mass was neglected. Preserving exact dependence on the muon mass in the case of muonic hydrogen 
we can present the proton polarizability contribution to the HFS in the form: 

AE P FS = - Zo T\ E F (A! + A 2 ) = (8* + 5 2 P )E F = <5 P E F , (8) 
27rm2(l + k) 

Ai=^°°^{^F 2 (Q 2 )/? (a)-4m3 jf° ^-(3 x {a, 9)G 1 (v, Q 2 )} , (9) 



™dQ 2 ' x 

where Vth determines the pion-nucleon threshold 



A 2 = -12m 2 / ^ / d^ 2 (a,0)G 2 ( l /,Q 2 ), (10) 
' 'o Q 2 Ju th 
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m 2 +Q 2 
2m 2 



(11) 



and the functions /?o,i,2 have the form: 



/3 (ct) = 2 , (7 = — 



(12) 



/MM) 



(t(-1 + ct6 



-2yT+CT + 2 - 4cr 2 - 2er 



VI + ct 



(13) 



+ 



2g(-6i 2 + e 3 / 2 VTT0 + + 2) 

y/0(l + 0) 



Q 2 " 



/MM) 



a(-l + a6) 



i - VT 



y/0{l + O) 



(14) 



IMQ 2 ) is the Pauli form factor of the proton and the proton anomalous magnetic moment k=1. 79284 7386(63) JTBfl . For 
past years there were not enough experimental data and theoretical information about proton spin-dependent structure 
functions. So, the previous study of the contribution AE^J Fg contains only estimation of the proton polarizability 
effects: 5 P ~ 1 + 2 ppm, or the calculation of main resonance contributions p^- |l4|Jl6] ] . The theoretical bound for 
the proton polarizability contribution to the HFS of electronic hydrogen is \S P \ < 4 ppm fl7fl . New estimation of the 
contribution S p in the hydrogen atom was done in fig) ] on the basis of modern experimental data and theoretical results 
on the structure functions Gi^fa, Q 2 ). By virtue of the fact that the muon to electron mass ratio m M /m e =206. 768266 
the lepton mass dependent terms in the functions f3i give appreciable contribution to the S p . It is our purpose here to 
calculate the correction S p for muonic hydrogen with the account muon mass dependent terms in (8)-(10). Previously 
we have considered also the other possible source of small uncertainty in the HFS interval, connected with the hadronic 
vacuum polarization jl9| . 

The polarized structure functions gi(^, Q 2 ) and g2(z y , Q 2 ) enter in the antisymmetric part of hadronic tensor W^, 



describing lepton- nucleon deep inelastic scattering (DIS) 20 



W [sl 

¥¥ (IV 



¥¥ (IV 1 



(15) 



w|S=f- 



q^ 



Wi(i/, 



P • q 



P • q 
— 



W 2 (^,Q 2 ) 



(16) 



W [A] 



^gi{v, Q 2 ) 



[(P-q)S' 3 -(S-q)P^: 



; 2 (^Q 2 ) 

(P • q) 2 



(17) 



is the totally antisymmetric tensor in four dimensions, g\{v, Q 2 ) = TE^yG\(y, Q 2 ), g2(z y , Q 2 ) = ni2i' 2 G2(^, Q 2 ), 
P is the four-momentum of the nucleon, x = Q 2 /2m2^ is the Bjorken variable, S is the proton spin four-vector, 
normalized to S 2 = —1, q 2 = — Q 2 is the square of the four-momentum transfer. The invariant quantity P-q is related 
to the energy transfer v in the proton rest frame: P • q = m.^v. The invariant mass of the electroproduced hadronic 
system W is then W 2 = m 2 , + 2m<xV — Q 2 = m 2 , + Q 2 (l/x — 1). Here Wi and W2 are the structure functions for 
unpolarized scattering. In the DIS regime the invariant mass W must be greater than any resonance in the nucleon. 
The threshold between the resonance region and the deep-inelastic region is not well defined, but it is usually taken 
to be at about W 2 = 4 GeV 2 . 
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(a) 



(b) 



FIG. 1. Feynman diagrams for proton polarizability correction to the hydrogen HFS. 

Hadronic tensor W^j, is proportional to the imaginary part of the off-shell Compton amplitude for the forward 
scattering of virtual photons on nucleons: 7*N — > 7*N. The photon-nucleon interaction depends on the photon 
polarization as well as on the nuclcon one. This gives four independent hclicity amplitudes of the form M a b ;cc j, with 
a, b, c, d values for the helicities of the photon and nucleon initial and final states: 

Ml,l/2;L 1/2, Mx,_i/2;l,-l/2) M ,l/2;0,l/2) ^-X,X/2;0,-X/2' 

These components correspond to the four structure functions Wi, W2, gi, g2- All other possible combinations of initial 
and final photon and nucleon helicities are related to the above by time reversal and parity transformation. 

The proton spin structure functions can be measured in the inelastic scattering of polarized electrons on polarized 
protons. Recent improvements in polarized lepton beams and nucleon targets have made it possible to make accurate 
measurements of nucleon polarized structure functions gi j2 in experiments at SLAC, CERN and DESY ||2l|-P7|l . The 
spin dependent structure functions can be expressed in terms of virtual photon-absorption cross sections ppf: 



l(^,Q 2 ) = 



m 2 • K 



87r 2 a(l + Q 2 /^ 2 ) 



CTl/2(^Q 2 ) -CT3/2(^,Q 2 ) + 



cttl(^, Q 2 ) 



(18) 



m 2 • K 



87r 2 a(l + Q 2 /i/ 2 ) 



-0-1/2 (f, Q 2 ) +03/2(2/,' 



-^=cttl(^,Q 2 ) 
v 4 



(19) 



where K = v — is the Hand kinematical flux factor for virtual photons, <J\/2i 03/2 are the virtual photoabsorption 
transverse cross sections for the total photon-nucleon helicity of 1/2 and 3/2 respectively, cttl is the interference term 
between the transverse and longitudinal photon-nucleon amplitudes. In this work we calculate contribution AE^ FS 
on the basis of the latest experimental data on structure functions gi, 2 (iA Q 2 ) and theoretical predictions for the cross 
sections 01/2,3/2/1% ■ 

III. POLARIZED STRUCTURE FUNCTIONS 
A. Deep-Inelastic region 



Our calculation of the contribution AE^ FS in the DIS region (W 2 > 4 GeV 2 ) is based as on the recent experimental 
data J2l|- |2^ , |2^| , so on the evolution equations for the polarized parton densities. The structure function gj is related 
to the polarized quark and gluon distributions by p8|,p9| 



l(x,Q 2 ) 



< e^ > 



[C N s ® Aq NS + C s ® AE + 2n f C g ® Ag] 



(20) 
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where < e 2 >= n f 1 Ya=i e ?' ® denotes convolution with respect to x. The nonsinglet and singlet quark distributions 
are defined as: 



i=l 



AC 1NS - E ( ~ 1 ) (Aqi + ' (21) 



where Aqi(x, Q 2 ) = q+(x, Q 2 ) — q-(x, Q 2 ) measures to what degree the parton of flavour q "remembers" its parent 
proton polarization, Ag(x, Q 2 ) is the longitudinally polarized gluon density, probed at a scale Q 2 . The evolution 
equations for the polarized parton densities are given by |p8| , p9| 

A AqNS = ^ P N| g) Aq NS , (22) 



d /A£\ 
dt \Ag) 



q»(t) 

2tt 



po po 

ss as 



/AS 



where t = ln(Q 2 /AQ). The coefficient functions C and the polarized splitting functions P are now known at NLO. 
To mitigate possible higher twist contributions we used the Fortran program for solving the Q 2 evolution equations, 
suggested in [^9|, only in the region of Q 2 > 1 Gev 2 . The comparison of obtained results for the polarized structure 
function gi(x, Q 2 ) with experimental data pl|- p3| , p6[ is presented in Fig. 2-4 at different Q 2 . Recent experimental 
data pl|- p3y2^ ] show that there are still large experimental errors for the polarized structure function gi(x, Q 2 ). So, 
application of the evolution equations in the nonresonance region of Q 2 > 1 GeV 2 allows to decrease substantially the 
theoretical errors in the 8 P calculation. In the region of Q 2 < 1 Gev 2 we used the experimental data for gi(x, Q 2 ). 
All of the data, including the SMC data at Q 2 < 1 GeV 2 , were fitted by the parameterization: 



gi(x, Q 2 ) = aix a2 (l + a 3 x + a 4 x 2 )[l + a 5 f(Q 2 )]F!(x, Q 2 



(23) 



where Fi = Wim2. The coefficients of the fits and different models for the form of the Q 2 dependence may be found in 
pjj . Numerical integration in (8) was performed with the f(Q 2 ) = — InQ 2 (fit IV), corresponding to the perturbative 
QCD behaviour. We have extrapolated relation (23) to the region near Q 2 = 0. Calculation of the second part of the 
correction <5 P in (8) for the nonresonance region was performed by means of the Wandzura-Wilczek relation between 
spin structure functions gi(x, Q 2 ) and g2(x, Q 2 ): 



^ w (x,Q*) = - gl (x,Cm 



L (t,Q 2 



dt 



,ww 

=2 



(24) 



Higher twist terms contribute to g2(x, Q 2 ) as well, but they are small enough. So when the transverse asymmetry Aj 
is not measured or is poorly known, g2 = is often used. 
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FIG. 2. Plots of gi(x,Q 2 )/Fi(x, Q 2 ) for Q 2 = 1 GeV 2 . The solid and dashed curves are the solution of the DGLAP evolution 
equation and the experimental fit correspondingly. Experimental points with total errors are taken from the paper [21]. 
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FIG. 3. Proton structure function xgi(x, Q 2 ) for Q 2 = 2 GeV 2 in the nonresonance region (solid curve). The dashed curve 
is the experimental fit. Experimental points correspond to the paper [21] 
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FIG. 4. Proton structure function xgi(x, Q 2 ) for Q 2 = 5 GeV 2 in the nonresonance region (solid curve). The dashed curve 
is the experimental fit. Experimental points correspond to the paper [21]. 



B. Resonance region 

The proton polarizability contribution to HFS in the resonance region is determined by the processes of photo- and 
electroproduction on nucleons of the pions and some prominent baryon resonances. The amplitudes of such reactions 
are shown on Fig. 5. 

To obtain correction (8) at the resonance region (W 2 < 4GeV 2 ) we use the Breit-Wigner parameterization for 
the photoabsorption cross sections in (18), (19), suggested in |^fj| -|36f. There are many baryon resonances that 
give contribution to photon absorption cross sections. We take into account only five most important resonances: 
P 33 (1232), Sn(1535), Di 3 (1520), Pn(1440), F 15 (1680). Considering the one-pion decay channel of the resonances, 
the absorption cross sections <J\/2 and may be written as follows p^j3^ ]: 

_ /k R \ 2 w 2 r 7 r R ,^ 4m p 2 

^1/2,3/2 -[^j (W 2 _ M 2 )2 + W 2 r 2 ot MrTr l A l/2.3/ 2 | (25) 

where Aj/2,3/2 are transverse electromagnetic helicity amplitudes, 

/ k V 1 /k 2 +X 2 \ j2 
r 7 =r R (-j (^ TW ) , X = 0.3 GeV. (26) 
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FIG. 5. Feynman amplitudes for proton polarizability correction in the resonance region. Solid, double solid, wave and 
dashed lines correspond to nucleon, baryon resonance, photon and pion correspondingly. 



The resonance parameters Tr, Mr, ji, ja, 
parameterization of one-pion decay width is 



rtot were taken from 15 3q|. In accordance with Refs. 



the 



■ R— »N 



M 



\QB. 



c 5 



c 2 



for the P 33 (1232) and 



R— »N: 



•(q) 



21+1 



q| 



C = 0.3 GeV 



i+i 



S 2 



(27) 



(28) 



for D 13 (1520), Pn(1440), F 15 (1680). 1 is the pion angular momentum and S 2 = (Mr- m p - m,) 2 + r R /4. Here q 
(k) and qR (kR) denote the c.m.s. pion (photon) momenta of resonances with mass M and Mr respectively. In the 
case of 5n(1535) we take into account 7rN and 77N decay modes |34| , [38f : 



r 



R— *7T,?7 



q 



q 2 

Drr^t R — ~ 



C 



7T ,7] 



p2 

7T,?7 



(29) 



where b,^ is the 7r (rj) branching ratio. 

The cross section cttl is determined by an expression similar to (25), containing product (Sy 2 ■ Ajya + A^Sjya) 
pl[ . The calculation of helicity amplitudes Ai/ 2 , A 3 / 2 and longitudinal amplitude S1/2, as functions of Q 2 , was done 
on the basis of constituent quark model (CQM) in |39|-|44|]. In the real photon limit Q 2 = we take corresponding 
resonance amplitudes from [g5|. At low center of mass energies the excitation of A(1232) resonance is of particular 
importance. At small Q 2 it is dominated by a magnetic dipole transition. The Q 2 dependence of A!/ 2 (Q 2 ), A 3 / 2 (Q 2 ) 
was extracted from the analysis made in Helicity amplitudes of the other resonances were taken from the 

calculations in CQM pl|-p4|. We have considered Roper resonance Pn(1440) as an ordinary qqq state fig} ]. As it 
follows from predictions of the quark model, the helicity amplitudes, which may be suppressed at Q 2 = 0, become 
dominant very rapidly with increasing Q 2 . 

The two-pion decay modes of the higher nucleon resonances (Sn(1535), Di 3 (1520), Pn(1440) and Fis(1680) were 
described phenomenologically using two-step process as in fl32| ]. The high lying nucleon resonance R can decay first 
into N* (P 33 (1232) or Pn(1440)) and a pion or into a nucleon and p- or er-meson. Then the new resonances decay 
into a nucleon and a pion or two pions: 



R 




7T — > N + 7T + 7T, 
N -> N + 7T + 7T. 



The total decay width of such processes can be presented as a phase space weight integral over the mass distribution 
of the intermediate resonance r= N*, p, a (a — n, N): 



- R— >r+a 



(W) 



P2tt 



W-m a 



dp • Pf 



M 2r r,tot(A«) 



(Mr - m a - 2m,) 2 + C 2 



7T (m 2 - m 2 ) 2 + M 2 r 2 tot (M) (W - m a - 2m,) 2 + C 2 



(30) 
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where C = 0.3 GeV, the factor must be taken from the constraint condition: rR^ r+a (Wrt) coincides with the 
experimental data, pf is the three momentum of resonance r in the rest frame of R. r r to t is the total width of the 
resonance r. The decay width of the meson resonance is parameterized similarly to that of the P33(1232): 

r(,)=r r ^(A) 2J ' +1 4±|!, g = 0.3 GeV, (31) 

where m r and fj, are the mean mass and the actual mass of the meson resonance, q and q r are the pion three momenta 
in the rest frame of the resonance with masses [i and m r , J r and T r are the spin and decay width of the resonance 
with mass m r . 

To calculate the nonresonant contributions to cross sections CT1/2, c 3 / 2 , 0"tl in the resonance region we have used 
predictions of the unitary isobar model (UIM) p6| , fi7| , [48] . Nonresonant background contribution to gi,2( x , Q 2 ) f° r 
W 2 < 4 GeV 2 can be derived from effective Lagrangians describing the electromagnetic 7NN, 77T7T interactions and 
the hadronic 7rNN system These Lagrangians lead to the well-known expressions for the CGLN amplitudes 

Fi, . . . ,Fe which give the dominant part of the background. The other part is related to vector meson exchange 
contributions. Within UIM taking into account Born terms, vector meson terms and the interference resonance 
background contributions to the polarized structure functions we can calculate the single-pion production am plitudes 



(see Fig. 5) using the on-line version of the numerical program MAID: http: / / www.kph-uni-mainz.de /MAID . We used 
it to consider nonresonance pion electroproduction contributions to Ai j2 . 

The Breit - Wigner five resonance parameterization of photon cross sections and constituent quark model results 
give good description of proton polarized structure functions in the resonance region. But the still existing difference 
between this description of %\,i{y, Q 2 ) and experimental data Jl^,[37j requires further improvement in the construction 
of spin dependent structure functions. It can be done considering contributions of additional baryon resonances in 
the large W range: S 3 i(1620), F 37 (1950), D 33 (1700), P 13 (1720), F 35 (1905) and accounting for different decay modes 
of such states. The particular significance in the study of the spin-dependent properties of baryon resonances belongs 
to Gerasimov-Drell-Hearn (GDH) sum rule |Q 

k 2 1 r°° du 

~T^ = T2- -K 2 M)-<t 3/2 M)] (32) 
4m 2 : 8n 2 a J Uth v 

The GDH sum rule rests on the basic physical principles and an unsubtracted dispersion relation applied to the 
forward Compton amplitude. Recently there was obtained first experimental data for the contribution to GDH sum 
rule in the energy range 200 < E 7 < 800 MeV [|TJ. Our expressions for the polarized structure functions gi,2(x, Q 2 ) 
show that sum rule (32) is valid with high accuracy. Here we must call special attention to the fact that the most 
important region for the integrals (9) and (10) on Q 2 variable lies from 0.1 Gev 2 to 1 GeV 2 . 

The second part of (9) gives especially large negative contribution to the correction Sf in the range of small Q 2 , 
where the contribution of A isobar is dominant. With increasing Q 2 its value falls and the total correction Sf has 
positive sign. 



IV. NUMERICAL RESULTS 



The processes of electromagnetic production of light quark baryon resonances represent the unique sphere for the 
investigation of nonperturbative properties of quantum chromodynamics (QCD). The direct test of the hadron quark 
model, chiral perturbation theory (ChPT), UIM, QCD sum rules and other approaches is connected with the study 
of the GDH sum rule pH , helicity amplitudes A1/2, A 3 / 2 , S;iy 2 and transition form factors [^2|,^3); the nucleon spin- 
dependent structure functions |2l|-p3|. The HFS measurement in the hydrogen atom and in the muonic hydrogen 
can be considered also for the verification of different effective field theories, because AEnFs(ep) is one of the most 
precisely measured quantities in atomic physics. 

In this work for the ground state hyperfine splitting interval we considered the part of the strong interaction effects, 
which manifest itself in the proton polarizability correction. The calculation of 5 P was based on three main ingredients: 

1. The evolution equations for the spin dependent structure functions in the nonresonance region, 

2. The phenomenological hadron quark model, isobar model in the resonance region. 

3. The experimental data for the nucleon polarized structure functions obtained at SLAC, DESY, CERN. 

The values of contributions <5f , 6% and the total contribution S p , obtained after the numerical integration in the 
resonance and nonresonance regions are as follows: 

Sl rcs = 3.13 • KT 4 , jf inonrcs = 2.01 ■ KT 4 , 5f = 5.14 ■ KT 4 , (33) 
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= -0.56-10 



-4 



2,nonres 



0.02 • 10~ 4 , 5% 



0.58- 10 



4 



(34) 




= 4.6 ±0.8 • 10 



(35) 



where the error, indicated in the expression (35), is determined by three main factors, connected with the polarized 
structure functions. We solved DGLAP equations in the NLO approximation, so possible uncertainty in 5 P can 
comprise near 10% of obtained result. The other source of the theoretical uncertainty arises from the experimental 
data errors in the Q 2 < 1 GeV 2 region. We estimated it at a level of about 20 % of the contribution i5 p at Q 2 < 1 GeV 2 
in the nonresonance region. Finally, the third part of uncertainty in (35) is controlled by the predictions of the 
constituent quark model for the electroproduction amplitudes A^, A 3 / 2 , Si/ 2 ||- We supposed that this one can 
reach 20 % of the contribution 5 P in the resonance region, comparing expressions (18)-(19) with the experimental 
data in the resonance region Jig] . On our opinion this is the main source of theoretical uncertainty now. The first way 
to obtain more reliable result in the resonance region can be based on using transition form factors of the nucleon to 
the baryon resonances with different values of J p , calculated on the basis of QCD sum rules as in |l| and of ChPT 
approach. The second one is connected with the progress in measurement of polarized structure functions gi j2 (x, Q 2 ) 
and perhaps lattice calculations. 

Muon mass dependent terms in relations (8)-(10), which are negligible for electronic hydrogen, give contribution to 
S p at a level of 25 %. In the case of muonic hydrogen E F (/xp) = 182.443 meV and the value of proton polarizability 
contribution to the HFS is equal to 0.084 meV for n=l and 0.011 meV for n=2 state. Last value can be important for 
the determination of the Lamb shift from the 2P-2S frequency measurement JjJ. The total correction (35) must be 
taken into account if the measurement of the ground state HFS in muonic hydrogen with the accuracy 10 became 
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